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ABSTRACT 

Several different processes could be changing the density in the core of a 
neutron star, leading to a departure from (3 equilibrium, quantified by the chem- 
ical potential difference = /i„ — /ip — /ig. The evolution of this quantity is 
coupled to that of the star's interior temperature T by two functions that quan- 
tify the rate at which neutrino-emitting reactions proceed: the net reaction rate 
(difference between (3 decay and capture rates), Fnetl^, ^/^), and the total emis- 
sivity (total energy emission rate in the form of neutrinos and antineutrinos) , 
etot(7', Here, we present a simple and general relationship between these 
variables, detot/ 96 fi = SFnet, and show that it holds even in the case of super- 
fluid nucleons. This relation may simplify the numerical calculation of these 
quantities, including superfluid reduction factors. 

Subject headings: dense matter — neutrinos — stars: neutron 



INTRODUCTION 



The simplest weak interaction process that could proceed in the core of a neutron star is 
the so-called direct Urea (or Durcd) process, which consists of the two successive reactions, 
/3 decay and capture. 

n ^ p + e + Ve p + e ^ n + Ue- (1) 

It is the most powerful of the neutrino processes potentially leading to the cooling of the 
neutron star. However, these reactions are allowed only if each of the Fermi momenta of 
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neutrons (n), protons (p), and electrons (e) is smaller than the sum of the two others (triangle 
condition). This implies a (highly uncertain) threshold in the matter density for the direct 
Urea processes. Another reaction, which overcomes this restriction, is the so-called modified 
Urea (or Mured) process, which involves an additional spectator nucleon: 

n + N^p + e + i>e + N p + e + N^n + u^ + N, (2) 

where the additional nucleon N can be either a neutron or a proton. The matter is trans- 
parent to neutrinos, which escape freely, transporting their energy away from the star. 

Furthermore, these two reactions bring nucleons into the state of beta (or chemical) 
equilibrium, which determines the concentration of neutrons and protons. The beta equilib- 
rium condition /i^ — /J'p + A*e involving the chemical potentials of the constituent particles 
implies the equality of Pe and Ue emission rates (or the net reaction rate set equal to zero). If 
a slight departure from equilibrium S/i = fj,n — [Jp — /ie is produced by any external or macro- 
scopic phenomenon that changes the density of a fluid element, one of the two reactions in 
equations (1) or (2) becomes more intense and changes the fraction of protons and neutrons 
towards the new equilibrium values (Le Chatelier's principle). In the subsequent evolution, 
the chemical imbalance affects the stellar interior temperature T by increasing the phase 
space available to the products of the neutrino-emitting reactions and by converting chemi- 
cal energy into thermal energy, and the temperature affects the chemical imbalance by also 
determining the rate at which reactions proceed. Several authors have investigated external 
processes that induce non-equilibrium /3 reactions, such as radial pulsation (Finzi 1965; Finzi 
& Wolf 1968), gravitational collapse (Haensel 1992; Gourgoulhon & Haensel 1993), a chang- 
ing rotation rate (Reisenegger 1995, 1997; Fernandez & Reisenegger 2005; Reisenegger et al. 
2006), and a hypothetical time- variation of the gravitational constant (Jofre et al. 2006). 
However, none of these has considered an adequate model for the effects of the likely Cooper 
pairing (superfluidity) of nucleons (see Reisenegger 1997 for a very rough estimate of their 
likely importance). 

Two basic rates are relevant in order to follow the coupled evolution of T and 6fi (see, 
e. g., Reisenegger 1995; Fernandez & Reisenegger 2005): the total (neutrino and antineutrino) 
emissivity (energy per unit volume per unit time), etot = + ^u, and the net reaction rate 
(number of reactions or emitted lepton number per unit volume per unit time), Fnet = Fp— F,^. 
Recently, Villain & Haensel (2005) analyzed the effect of nucleon superfluidity on the net 
reaction rates, Fnet, calculating superfluid reduction factors for the direct and modified Urea 
processes by means of sophisticated numerical methods. Nevertheless, they did not evaluate 
the total emissivity of neutrinos, etot, which is also required in order to compute the eventual 
time-evolution of T and S/i. 

In this work, we present a curious relationship between etot and Fnet, both considered 
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as functions of the chemical imbalance parameter S/j,, that holds even in the superfluid case: 

The advantage of this relation is that, having computed one of the functions numerically 
(as done by Villain & Haensel 2005), the other is obtained very easily, saving time in the 
calculation. 

In the non-superfiuid case, Reisenegger (1995) calculated etot and Fnet analytically. For 
Durca processes. 



^Durca / 
-tot 



1071^2 + 315m^ + 21 



(T, 5/.) = 6? ( 1 + "::r "^^^ ) , (4) 



457 



and 



.Du.ca,^ ,D714ii^ + 420x.^ + 42t.« 



r^r^'iT, 6^1)6^1 = , (5) 

where u = 5ijl/{'kT) is a dimensionless parameter, and = eJ^"'^'^^(T', 0) oc is the equilib- 
rium emissivity. For Murca processes, 

,--(T. S,) = 4- (l + " + ^74^^^°"° + ""' ) . (6) 

and 

iioio 

with = e^"'''^'^(T, 0) oc (see Haensel 1992 for precise estimates). In these two pairs of 
polynomial expressions, it is easy to verify our proposed relation (eq. 3). The purpose of this 
paper is to show that it is valid beyond these simple cases, encompassing, for example, the 
cases of superfluid neutrons and/or protons. A brief and clear discussion of superfluidity in 
neutron star was given by Yakovlev (2001), and a rigorous and deep analysis can be found 
in Lombardo & Schulze (2001). 



2. DERIVATION 

We now show that equation (3) holds regardless of the nucleons being superfluid or not. 
Along this paper we will use natural units, with H—kB — c—l. 



We express the phase space factors as 

dpi = pldpi dQi = D{Ei)<iEi dQj, 



(8) 
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where Ei{pi) is the particle energy,^ and 

(9) 

is the density of states and dilj is the sohd angle element in the direction of pi. In order to 
show the generality of our derivation, we left the electron and nucleon densities of states^ 
expressed implicitly in the following calculation. For neutrinos, we will need the explicit 
expression 

dp,^ EldE^dQ,, (10) 
and we may assume neutrino isotropy, so 

dp^^AnEldE,. (11) 



In the interior of the neutron star, the temperature is much smaller than the Fermi tem- 
perature. On the other hand, the neutrino momentum is proportional to the temperature of 
the star, and the other momenta are essentially their respective Fermi momenta (degenerate 
matter). Therefore, we assume that \pi,\ <^ |pi| {i = n, p, e), and neglect in the Dirac 
delta function of momentum. This approach will allow us to integrate more easily over the 
orientation of the neutrino momentum. 

Using the following dimensionless variables: 

El- Pi ,. X Ey 

Xi = — - — [i = n, p, e), Xy = — and u= — , (12) 
the total emissivity etot can be expressed as 

6tot(T,<5/.) = ^(|M|^)n/_ (13) 

(Yakovlev et al. 2001, eq. 115), where (|Mp) is the squared transition amplitude for the 
f3 decay and capture processes, averaged over initial and summed over final spin states; 
furthermore averaged over the direction of the neutrino momentum, this results in a roughly 



•^The dispersion relation for the superfluid nucleons can be written as E{p) = jj + sgn(y) — 
Pf)\/v'p{p — PfY + where pp denotes the Fermi momenta. The Cooper pairing occurring around the 
Fermi surfaces, allows us to approximate (l/2m)(p^ ~ Pf) — v\{p — pp^ , where vp = pp/m* and nf are 
the Fermi velocity and effective particle mass for superfluid nucleons, respectively. 

"^The density of states for superfluid nucleons is D{E) = ^[^"^'j^^^, Q{\E - fj,\ - A). 
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constant value, which can be taken out of the integrals (for more details see Shapiro & 
Teukolsky 1983 or Yakovlev et al. 2001): 

(|Mp) ~2G^cos2^c(l + 35i), (14) 

where Gp, 9c, and qa are the Fermi weak interaction constant, the Cabibbo angle {sm9c = 
0.231), and the Gamow- Teller axial vector coupling constant, respectively. The operator fl 
contains the integrals over the orientations of the particle momenta 



fi = 47r JJJ dilnd^pd^eS{Pn-Pp-Pe): (15) 

and /_ integral includes the integrations over dimensionless particle energies: 

/oo poo roo POO 

dXeDefe dXnDnfn dXpDpfp / dx^X^ 
-oo J —oo J —oo Jo 

X [5{Xn + Xp + Xe — Xi,+ u) — S{Xn+Xp + Xe—X„ — u)] , (16) 

where the /j's (for i — p, n, e) are the Fermi-Dirac distributions, /j = (ef + l)""*^, for 
some of whose arguments the signs have been redefined when using the symmetry property 
fi{x) = 1 — fi{—x). The functions D{xi) are approximately symmetric around x = for all 
relevant cases. 

When considering situations with neutron or proton superfluidity, the possible anisotropy 
of gaps (appearing in the density of states function) in the expression 7_ should be taken 
into account in the preceding Q integration. 

Similarly, the net reaction rate is 

^5 



V^,,{T,5i,) = -—A\M\')ni+, (17) 



(27r) 



where 



X) POO poo poo 

dXeD{Xe)fe dXnD{Xn)fn dXpD{Xp)fp / dx^X^ 
oo J —oo J —oo Jo 

X [6{Xn + Xp + Xe — Xt,+ u) + 5{Xn + Xp+Xe — X„—u)] . (18) 

Note that the integral now contains only two powers of Xj^. 

The total emissivity (given by eq. [13]) depends on the chemical imbalance Sn only 
through the dimensionless non-equilibrium parameter u contained in in the argument of 
the Dirac delta function. In order to calculate the derivative deiot/dSn, we define 

Z± = Xn + Xp + Xe — Xi,± U, (19) 
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so we can rewrite the derivative 

^ [6{z^)-6{z.)] = A [-6{z^)-6{z.)] (20) 

and do an integral by parts in order to obtain the derivative of equation (16)^: 

-g^= dXeD{xe)fe dx„D(x„)/„ / dXpD{xp)fp dx^,{3xl) [S{z+) + 5{z^)] ^ 31+. 

(21) 

Thus, in a straightforward way, we reach the proposed relation (eq. 3): 

It is straightforward to verify that an analogous derivation can be made for modified 
Urea processes. 



3. COOLING VS. HEATING 

As a useful application of the proposed relation, we analyze the balance between heating 
and coohng due to /3 processes. 

The net local heating rate can be written as 

^net — Tnet*^/^ ~ ^tot; (22) 

where the first term is the rate of dissipation of chemical energy, and the second is the energy 
loss rate due to neutrino emission. Using our relation (eq. 3), which links etot and Fnet, we 
obtain 

Sl^^ d / etot 
3 86 /J \Sfi^ 

If etot increases faster than the net heating is positive (/inet > 0), and viccvcrsa. 



The net heating rate can also be written as 



''«. = %ifl^-3). (24) 



3 \dhiSfj, 

Without Cooper pairing, and in the limiting case oi S/i ^ T, we know that etot oc S/i^ for 
Durca processes, and etot oc 5//^ for Murca processes. For these cases, we easily reobtain that 



^The boundary term is of the form x^/{l + e^)|o° and vanishes in both Umits. 
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the total energy released, Fnetf^A*, is distributed in fixed fractions among internal heating, /inet, 
and neutrino emission, as /inet = ^tot = ^^net^lJ' for Durca, and /inet — |rnet<^A*) ^tot — |rnet<^A* 
for Murca processes (Fernandez & Reisenegger 2005). 

4. CONCLUSIONS 

We have proven a simple, general relationship (eq. 3) between the main rates character- 
izing non-equilibrium (3 processes in both superfiuid and non-superfluid neutron star matter. 
This relation could simplify the evaluation of these quantities in superfiuid neutron star 
models, complementing numerical calculations such as those of Villain & Haensel (2005). 
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